A unified theory of the Elliott- Yafet and the D'yakonov-Perel' spin-relaxation mechanisms 



(N 

o 

(N 

> 
O 

in 



Peter Boross, 1 ' 2 Balazs Dora, 1 ' 3 Annamaria Kiss, 1 ' 4 and Ferenc Simon 1 ' 5 '© 

'Department of Physics, Budapest University of Technology and Economics, Budapest, Hungary 
2 Department of Materials Physics, Eotvos University, Budapest, Hungary 
3 BME-MTA Exotic Quantum Phases Research Group, Budapest University of Technology and Economics, Budapest, Hungary 
4 Wigner Research Centre for Physics of the Hungarian Academy of Sciences, Budapest, Hungary 
5 Condensed Matter Physics Research Group of the Hungarian Academy of Sciences, Budapest, Hungary 

(Dated: November 6, 2012) 

We present a unified treatment of the Elliott- Yafet (EY) and the D'yakonov-Perel' (DP) spin-relaxation mech- 
anisms using the Mori-Kawasaki formula, which gives the spin-relaxation rate to lowest order in the spin-orbit 
coupling (SOC) but to infinite order in the quasi-particle scattering rate, V. We consider a four-state Hamiltonian 
of the conduction and a nearby (valence) band with spin degeneracy, including SOC between adjacent bands 
(inter-SOC) and within the same band (intra-SOC). We find in agreement with the expectations that intra-SOC 
yields the DP- whereas the inter-SOC the EY-like result. However, we identify parameter domains of Y and 
the band structure where a crossover occurs between the two types of spin-relaxation mechanisms. The result 
ultimately connects the EY and the DP spin-relaxation mechanisms into a unified description and it leads to a 
better understanding of spin-relaxation in strongly correlated systems and where band degeneracy plays a role 
such as e.g. in graphene. 
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PACS numbers: 76.30.Pk, 71.70.Ej, 75.76.+j 

Introduction. Spintronics is an emerging new paradigm 
which intends to perform calculations and store information 
using the spin-degree of freedom of electrons and it could 
eventually replace conventional electronics (U]. A future spin- 
tronic device would operate with a spin-polarized ensemble of 
electrons whose spin-state is manipulated in a transistor-like 
configuration and is read out with a spin-detector (or spin- 
valve). Clearly, the utility of spintronics relies on whether the 
spin-polarization of the electron ensemble can be maintained 
sufficiently long. The basic idea behind spintronics is that co- 
herence of a spin-ensemble persists longer than the coherence 
of electron momentum due to the relatively weaker coupling 
of the spin to the environment. The coupling is relativistic and 
has thus a relatively weak effect known as spin-orbit coupling 
(SOC). 

The time characterizing the decay of spin-polarization is 
the so-called spin-relaxation time (often also referred to as 
spin-lattice relaxation time), r s . It can be measured either 
using electron spin-resonance spectroscopy (ESR) J3l or in 
spin-transport experiments . Much as the theory and ex- 
periments of spin-relaxation measurements are developed, it 
remains an intensively studied field for novel materials; e.g. 
the value of r s is the matter of intensive theoretical studies 
JH - 11 ] and spin-transport experiments 1 12 - 15 ] in graphene at 
present. 

Conventionally, the theory of spin-relaxation in metals 
and semiconductors is discussed along disjoint avenues, the 
so-called Elliott- Yafet (EY) and the D'yakonov-Perel' (DP) 
mechanisms due to reasons described below. The EY theory 
1 HI 12 ] describes spin-relaxation in metals and semiconduc- 
tors with inversion symmetry. Therein, the SOC does not spli t 
the spin-up/down states (|y), |4-}) in the conduction band 11811 . 
however the presence of a near lying band weakly mixes these 
states while maintaining the energy degeneracy. The nomi- 
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FIG. 1. Schematic relaxation processes: a) the Elliott- Yafet mecha- 
nism, b) the D'yakonov-Perel' mechanism. Note that spin-scattering 
occurs rarely in the EY scenario (typically for every 10 4 ..10 6 th scat- 
tering in alkali metals), whereas the spin direction continuously pro- 
cesses around the internal magnetic field due to SOC in the DP sce- 



nally up state reads: f) = cik |t) + ^fe 14-} (here a^, bk are 
band structure dependent) and b-^ja^ = L/A, where L is 
the SOC matrix element between the adjacent bands and A 
is their separation. E.g. in alkali metals L/A ps 10 -2 ..1CP 3 
[Ref.ll7|]. Elliott showed using first order time-dependent per- 
turbation theory that an electron can flip its spin with proba- 
bility (L/A) 2 at a momentum scattering event. As a result, 
the spin scattering rate (F s = K/2t$) reads: 
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where V = H/2t is the quasi-particle scattering rate with r 
being the corresponding momentum scattering (or relaxation) 
time. This mechanism is schematically depicted in Fig.Hk. 

For semiconductors with zinc-blende crystal structure, such 
as e.g. GaAs, the lack of inversion symmetry results in an 
efficient relaxation mechanism, the D'yakonov-Perel' spin- 
relaxation Il9ll . Therein, the spin-up/down energy levels in 
the conduction bands are split. The splitting acts on the elec- 
trons as if an internal, fc-dependent magnetic field would be 
present, around which the electron spins precess with a Lar- 
mor frequency of fi(fe) = J£(k)/K. Here .Sf(fc) is the energy 
scale for the inversion symmetry breaking induced SOC. Were 
no momentum scattering present, the electron energies would 
acquire a distribution according to Kl(k). In the presence 
of momentum scattering which satisfies Sl(k) ■ t -C 1, the 
distribution is "motionally-narrowed" and the resulting spin- 
relaxation rate reads: 



r s 



s.DP 



r 



(2) 



This situation is depicted in Fig. QJ>. Clearly, the EY and DP 
mechanisms result in different dependence on V which is of- 
ten used for the empirical assignment of the relaxation mech- 
anism lEoll . 

The observation of an anomalous temperature dependence 
of the spin-relaxation time in MgB2 12 ill and the alkali ful- 
lerides 12211 and the development of a generalization of the 
EY theory highlighted that the spin-relaxation theory is not 
yet complete. In particular, the first order perturbation theory 
of Elliott breaks down when the quasi-particle scattering rate 
is not negligible compared to the other energy scales. One ex- 
pects similar surprises for the DP theory when the magnitude 
of e.g. the Zeeman energy is considered in comparison to the 
other relevant energy scales. 

Herein, we develop a general and robust theory of spin- 
relaxation in metals and semiconductors including SOC be- 
tween different bands and the same bands, provided the crys- 
tal symmetry allows for the latter. We employ the Mori- 
Kawasaki theory which considers the kinetic motion of the 
electrons under the perturbation of the SOC. We obtain a gen- 
eral result which contains both the EY and the DP mecha- 
nisms as limits when the quasi-particle scattering and the mag- 
netic field are small. Interesting links are recognized between 
the two mechanisms when these conditions are violated: the 
EY mechanism appears to the DP-like when V is large com- 
pared to A and the DP mechanism appears to be EY-like when 
the Zeeman energy is larger than T. Qualitative explanations 
are provided for these analytically observed behaviors. 

Results and Discussion. The minimal model of spin- 
relaxation is a four-state (two bands with spin) model Hamil- 
tonian for a two-dimensional electron gas (2DEG) in a mag- 




FIG. 2. Color online. The band structure of a 2DEG in a magnetic 
field. The effects of the weak SOC are not shown. Vertical arrows 
show the energy separations between the relevant bands. 



netic field, which reads: 



n — Ho + Hz + %scatt + rl 

%> = 6k ' a C k,a,s C k,a,s 

Hz = A z sc l a ,s c k, a , s 



so 



k,a.s 



H so = 



(3a) 
(3b) 

(3c) 

(3d) 



where a = 1 (nearby), 2 (conduction) is the band index with 
s = (t) , (I) spin, e fc , Q = h 2 k 2 /2m* a - <5 Q .iA is the single- 
particle dispersion with m* = (— l) a m* effective mass and 
A band gap, A z = g^B z is the Zeeman energy. T^scatt is re- 
sponsible for the finite quasi-particle lifetime due to impurity 
and electron-phonon scattering and i a;Q '. S;S ' (k) is the SOC. 

The corresponding band structure is depicted in Fig. [2] The 
eigenenergies and eigenstates without SOC are 



ek.a,s = £fe, Q + sA z (4a) 
|1,4)= [1,0,0, Of |l,t} = [0,1,0, 0] T (4b) 
|24) = [0,0,1,0] T |2,t) = [0,0,0,l] T . (4c) 

The most general expression of the SOC for the above lev- 
els reads: 
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where B± is the magnetic induction of the electromagnetic 
radiation, x± is the imaginary part of the spin-susceptibility, 
fXo is the permeability of vacuum, and V is the sample vol- 
ume. The spin-susceptibility is related to the retarded Green's 
function as 



TABLE I. Effect of the presence or absence of the inversion sym- 
metry on the intra- (Jzf) and inter-band (L) SOC and on the energy 
splitting of spin-states in the same band, et,^ — £k,i- 



(fc) 



/ £&ft J^?if X-f-f L±f \ 

"^^W- "^^X-X- "^''1*4' 

^TT ~^Tt 

V L n L u -&n ■Z'u J 



(5) 



-'a, a .s.s 



(fc) = 



(6) 



where ££ ss i (k), L ss i (k) are the wavevector dependent intra- 
and inter-band terms, respectively. The terms mixing the same 
spin direction can be ignored as they commute with the S z 
operator and do not cause spin-relaxation. The SOC terms 
contributing to spin-relaxation are 

/ JSf L \ 

JSf+ Lt 

L if 

V Lt J2ft o / 

Table U summarizes the role of the inversion symmetry on 
the SOC parameters. For a material with inversion symme- 
try, the Kramers theorem dictates (without magnetic field) that 

(fc) = ej. (k) and thus Jz? = 0, which term would otherwise 
split the spin degeneracy in the same band. When the inver- 
sion symmetry is broken, Jz? is finite and the previous degen- 
eracy is reduced to a weaker condition: (k) = ej_ (— k) as 
the time-reversal symmetry is retained. 

We consider the SOC as the smallest energy scale in our 
model {££ (hg) , L {hp)), while we allow for a competition of 
the other energy scales, namely Az, Y and A, which can be of 
the same order of magnitude, as opposed to the conventional 
EY or DP case. We are mainly interested in the regime of a 
weak SOC, moderate magnetic fields, high occupation, and a 
large band gap. We treat the quasi-particle scattering rate to 
infinite order thus large values of Y are possible. 

The energy spectrum of the spins (or the ESR line- 
width) can be calculated from the Mori-Kawasaki formula 
023112411 . which relies on the assumption that the line-shape is 
Lorentzian. This was originally proposed for localized spins 
(e.g. Heisenberg-type models) but it can be extended to itin- 
erant electrons. The standard (Faraday) ESR configuration 
measures the absorption of the electromagnetic wave polar- 
ized perpendicular to the static magnetic field. The ESR signal 
intensity is 



1(c) 



2flQ 



X 'i(q = 0,oj)V, 



(7) 



X 'i M = -ImG| +s _ H , (8) 

with S = S x ± iS y , from which the ESR spectrum can be 
obtained. 

The equation of motion of the S + operator reads as 



dS+ 
dt 



l -[H,S + ] = l -[U z ,S + ] + l -[Uso,S + l (9) 



i A 



where A = 4 \Hso : S + ] is the consequence of the SOC. The 
Green's function of S + S~ is obtained from the Green's func- 
tion of A* A as 



2(5 2 ) -([A(0),S-(0)]) + G% a (oj) 



ui 



Az 



Ai 

n 



r 



(10) 



The second term is zero without SOC thus a completely sharp 
resonance occurs at the Zeeman energy. The line-shape is 
Lorentzian for a weak SOC: 



G R S+S - M = 



2(S Z 



A z _ E(uO ' 



h h 

where the self-energy is 

-([A(0),S- (0)}) + G% a (u) 



E(w) 



2(S Z ) 



(ID 



(12) 



which is assumed to be a smooth function of u near Az /h. 
The spin-relaxation rate is equal to the imaginary part of 

£ (lj) as 



r s = 



hnG% A 



As 

, h 



(13) 



The G A t A (to) correlator is obtained from the Matsubara 
Green's function of A^A, given by 



Gat A ( il/ m) 



dre w ™ T (TtAHt) A(0)) . (14) 



The effect of ?4catt is taken into account in the Green's func- 
tion by a finite, constant momentum-scattering rate M25I1 . 

The most compact form of the spin-relaxation is obtained 
when the Fermi energy is not close to the bottom of the con- 
duction band (/i > A) and a calculation (detailed in the Sup- 
plementary Material) using Eq. ( fT3l leads to our main result: 



r s = 



4r|jsf (fcp)| 2 , 4r|L(fep)r 
4r 2 



A! 



4f 2 + A 2 (fc F ) ' 



(15) 
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FIG. 3. r s as a function of T according to Eq. J32| > concerning sep- 
arately the contributions due to intra- (upper) and inter-band (lower) 
SOC's. The insets are schematics of the band-structure including the 
broadening due to V (double arrows indicate the matrix elements). 
The arrows indicate the equivalence of the different spin-relaxation 
regimes. The conventional DP and EY scenarios are those in the 
top-right and bottom-left corners, respectively. 

Results in more general cases are discussed in the Supplemen- 
tary Material. 

The contributions from intra- (Jz? (fcp)) and inter-band 
(L (kp)) processes are additive to lowest order in the SOC 
and have a surprisingly similar form. A competition is ob- 
served between lifetime induced broadening (due to T) and 
the energy separation between states (A (kp) or Az). The sit- 
uation, together with schematics of the corresponding band- 
structures, is shown in Fig. [3] When the broadening is much 
smaller than the energy separation, the relaxation is EY-like, 
r s oc r, even when the intra-band SOC dominates, i.e. for 
a material with inversion symmetry breaking. This situation 
was also studied in Ref. |2g and it may be realized in III-V 

I 



semiconductors in high magnetic fields. For metals with in- 
version symmetry, this is the canonical EY regime. 

When the states are broadened beyond distinguishabil- 
ity (i.e. T ^ A (kp) or Az), spin-relaxation is caused 
by two quasi-degenerate states and the relaxation is of DP- 
type, r s oc 1 /r, even for a metal with inversion symmetry, 
Jz? = 0. The latter situation was observed experimentally in 
MgBoiRef.EIl) and the alkali fullerides (K 3 C 60 and Rb 3 C 60 ) 
(Ref. l22l) . which are strongly correlated metals with large T. 
When the intra-band SOC dominates, i.e. for a strong inver- 
sion symmetry breaking, this is the canonical DP regime. 

Similar behavior can be observed in other models (see 
Ref. I25I0 . r s oc r and T s oc l/T remain valid in the two 
different limits but the intermediate behavior is not universal. 

A particularly compelling situation is the case of graphene 
where a four-fold degeneracy is present at the Dirac -point and 
both inter- and intra-band SOC are present thus changing the 
chemical potential would allow to map the crossovers pre- 
dicted herein. 

Conclusions. In conclusion, we presented a unified treat- 
ment of the Elliott- Yafet the the D'yakonov-Perel' relaxation 
mechanisms. We found that depending on the quasi-particle 
scattering rate and the energy separation between states con- 
nected by the SOC terms, the spin-relaxation time crosses 
over smoothly from a DP to an EY-like behavior. In par- 
ticular, the r <C A (kp) region of the EY term corresponds 
to that of the DP in the T -C A^ region, since the lifetime 
broadening is small and the separate energy states can clearly 
be identified. As opposed to this, the T ^ A (kp) regions 
smears out the bands, connected by the SOC, causing them 
to be quasi-degenerate in terms of their energy. These ob- 
servations provide the ultimate link between these two spin- 
relaxation mechanisms, which are conventionally thought as 
being mutually exclusive. 
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SUPPLEMENTARY INFORMATION 



This Supplementary Material is organized as follows: we first discuss the technical details of the calculations, i.e. how Eq. (15) 
is obtained from Eq. (14). Second, we discuss the generalization of the spin-relaxation for two kinds of band dispersions 
(quadratic and linear or linearized) when the restrictions concerning the relative magnitude of the parameters (A,/z,m*) are 
lifted. We arrive at the overall conclusion that while the quantitative details of the V dependent spin-relaxation (r s ) are modified 
and no closed form of the result can be provided in the most general case, the overall trends, which characterize the EY and 
DP behaviors and in particular the crossover between the two, remain valid. Wherever possible, we provide closed form results 
though. Finally, we discuss the spin-relaxation for a model where only Rashba-type spin-relaxation is present. 
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I. TECHNICAL DETAILS OF THE CALCULATIONS 

We consider Eq. (14) of the main manuscript as a starting point.. The Matsubara Green's function of A' A can be written as 

Ga*A ( il/ m) = X] |(^ a > s\A^\k 7 a',s')\ 2 Q a:S (iu) n ,k)g a > )S < (ioj n +iu m ,k) , (16) 



where 



iu n ,k,a,a f ,s,s' 



Qa,s fan, k) 



\ (ek, a ,s - i") + ^ S g n 



(17) 



is the Matsubara Green's function of fermionic field operators in band (a) and spin (s). The effect of H sca tt is taken into account 
by the finite momentum-scattering rate, V. 

Using the relationship between the Green's function and spectral density, the Matsubara summation in Eq. ( fT6l ) yields 

I' 



1 2 /" /'OC 

Q A^ v m) = J- 2 X] |(fc,a,s|.4 t |fc,a',s')| / / efcj'cL; 

fe,a,a',s,s' J-ooJ-oo 



iv m — Ul" + OJ 



where 



p a . s (u>, k) 



-2 



A*)]' 



(19) 



is the spectral density. By taking the imaginary part after analytical continuation, the energy integrals can be calculated at zero 
temperature. Then, by replacing momentum summation with integration, we obtain 



A f°° 2 

lmG AfA {uj) = —J dkk ^ | a -> a '' s ') I £a,s, a >,s> (k, w) , 



(20) 



where 



4r 2 + (e fc , ajS - e k , a ', s ' - ^w) 



arctan 



arctan 



arctan 



arctan 



e klCe > s > + tuv 





"p 2 + r 2 






ek,a',s' - ftw) 


4r 2 + (e k , a , s - e k ,a',s' - ftw) 2 




(e k , a ,s - fw) + T 2 


(efc, Q ',s' + + T 2 



(21) 



and efe. QjS = ek.a.s — /i, A is the area of the 2DEG. The matrix elements of the A^ operator are 

(fc, a, s\A' I k, a', s') 



/'-if' —L \ 
L 

-L -if 
\ L if / 

We determine the expectation value of the z-component of electron spin following similar steps as 



(22) 



/CO (jjjj 1 
-co 



iuj n ,k,ck,s 



where 



(a,s (&) = 7: arctan 

Z 7T 



^ (A:,a,s|5 , 2 |A:,a,5)Ca,s (&) , (23) 



(24) 
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FIG. 4. Band structure of quadratic dispersion model. Vertical arrows show the energy separations between the relevant bands. 



The matrix elements of the S z operator are 



(k,a,s\S z \k,a,s) = - 



The spin-relaxation rate can be obtained as 



2{S Z ) 



( -1 \ 

10 

0-10 

V 1 J 



■p intra _j_ pinter 



(25) 



(26) 



We note this is the sum of intra- and inter-band terms which are described separately. 



II. SPIN-RELAXATION FOR DIFFERENT MODEL DISPERSIONS 



A. Quadratic dispersion model 

First, we discuss a quadratic model with the e^.a = f^k 2 /2m* a — <5 Q; i A single-particle dispersion. In the conduction band, 
the quasi-particles are electron-type (i.e. 7712 > 0) however the quasi-particles of a nearby band are hole-type (i.e. m\ < 0). 
The band structure is depicted in Fig.|U This model describes well two bands of the spectrum of semiconductors, however in a 
realistic case (e.g. for Si and GaAS) there are more nearby bands characterized by different band gaps and effective masses. 



1. The intra-band term 



An important and general limit of the model is when the Zeemann energy is much smaller than the band gap (i.e. Az <C A), 
and both the spin-up and spin-down states are occupied in the conduction band (i.e. Az <C /i). In this limit, the intra-band term 
can be expressed as 

pintra _ 

s " 4r 2 + Al • K n 

This term comes from processes within the conduction band and the nearby band does not give a contribution to the intra-band 
term. 
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2. The inter-band term in the general case 

In the Az -C /•«, A limit, the inter-band term can be determined but it takes a more complicated form. When the broadening 
is much smaller than the energy separation at the Fermi wavenumber (i.e. r« A (fcp) = A + (1 — m*,/m\) ff) the inter-band 
spin-relaxation has the form of 

Pinter _ 4 |Z, (fc F )| 2 

Fs " A 2 (fc F ) ' ( 8) 

which is directly proportional to the momentum-scattering rate. 

When the broadening is much larger than the energy separation (i.e. T ^> A (fc F )) the spin-relaxation reads 

winter -^m\m\ \L (fcp) | 2 1 



K - ml) 



2 



V 



(29) 



which is inversely proportional to the momentum-scattering rate. 



3. The inter-band term in the case ofm* = — m*t 
If the effective masses in the two bands have different signs but the same magnitude, the spin-relaxation rate is obtained as 



4|L(fc F )[ 2 r 
4r 2 + a 2 (fc F ) 



rin r2 +f 2 



(A + 2fi) ( 7r + arctan & - arctan 



(30) 



4. The inter-band term in the case ofm* = — m*\ and /! > A 

When the Fermi energy is not close to the bottom of the conduction band, the logarithmic term can be neglected and we obtain 
the most compact form of the inter-band spin-relaxation rate as 

into= 4|£(fc F )| 2 r 

4r 2 + A 2 (fc F )' K ' 

5. Summary of the result for the quadratic model 
Summing of intra- and inter-band term yields Eq. (15) of the paper: 

4|if(fc F )l 2 r 4|£(fc F )l 2 r 
4r 2 + A 2 4r 2 + A 2 (fc F )' K ' 

This is the main result, which is presented in the manuscript. 

We note that similar expressions can be obtained if the chemical potential lies in the nearby band. 

B. Linear band dispersion models 

Herein, we discuss the spin-relaxation time for linear model dispersions. The importance of studying this model is two-fold. 
First, every non-linear band dispersions can be linearized at the Fermi wave-vector and the plausible expectation is that the 
spin-relaxation rate can be obtained as a sum of the linearized segments. Second, spin-relaxation can be calculated for the linear 
band dispersion model and as we show below, the qualitative result, i.e. dependence of r s on T for the intra- and inter-band 
processes, is unchanged compared to the quadratic band dispersion even if the numerical factors are different. This proves that 
our calculation of the spin-relaxation is robust against the details of the band dispersion. 

The linear band-dispersion models can take two characteristically different scenarios: those with lines with slopes of the 
opposite and the same sign. The situation is depicted in Fig. [5] The first situation can be obtained e.g. from linearizing a 



quadratic band dispersion in Fig.|4]the Fermi wavenumber and the second occurs e.g. for MgB 2 as shown in Ref. 12111 . 
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1. Linear band dispersions with the opposite slope 

We consider that the higher lying conduction band has a positive Fermi velocity thus e'™ 2 = hv^ (k — kp) + \i, where the 
Zeeman-energy is neglected. The nearby valence band can be approximated with a flat band with zero Fermi velocity: e\\ = 
—A + [ira* 2 lm\, where the Zeeman splitting is also neglected too. 

Then, our calculation yields for the intra-band contribution to the spin-relaxation rate: 

pintra _ ^ \% Q*t m , 

s ~ 4r 2 + A 2 ' 

Similarly, we obtain for the inter-band term: 

pinter = 4|L(fc F )| 2 r 

T 2 + A 2 (kp) ' { } 

which look likes as if it was obtained from the quadratic model Eq. ( |3ll except the multiplication factor of the T 2 in the 
denominator. This is the result that we considered a zero Fermi velocity of the valence band. 



2. Linear hand dispersions with the same slope 

The second linear model has two linear bands (apart from the spin) with positive Fermi velocities of different magnitudes. 
The two bands cross the Fermi level at two separate points. The band structure of this model is depicted in Fig. [Sj?. This model 



describes the spectrum around the Fermi energy in e.g. MgB 2 as it was shown in Ref. |21 
The intra-band term is similar to the previous results and it reads: 



The inter-band term can be expressed as 



intra _ 4r|if(fc P )r 

4r 2 + A 2 ■ ; 



ter 4r|£(fc F )l 2 

X S = IA 4- A ^ 2 ' 
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where Ai and A2 are the distances of the two bands when one of the bands cross the Fermi level. The formula is symmetric in 
these two variables which means that the two bands change their roles as conduction and valence bands for the two Fermi level 
crossing points. 

A special case is when = uf, i.e. when the two linear bands are parallel therefore Ai = A2 = A. This yields a result 
which is similar to the case of the quadratic dispersion and reads: 

pinter = 4r|£(fc F )| 2 

4r 2 + A 2 ' K ' 



III. THE SPIN-RELAXATON FOR A MODEL WITH RASHBA-LIKE SOC 



Now we determine the spin-relaxation rate of a model where only Rashba-type SOC is present. The Rashba-like SOC can be 
written as 



Hso = fiA (a x k y - cr y k x ) , (38) 

/v.CV 

where a x , a y are the Pauli matrices. The matrix element of intra-band SOC can be expressed as Jz? (k) = HXk. We can expand 
Eq. (O to get 

8m$ M A 2 

S ~4r2 + Ai- 

Using the spin and momentum life-times instead of relaxation-rates, we obtain 

1 8m2^X 2 t 



h 2 1 , fAzl 



(40) 



A similar expression was obtained recently (Eq. (40) in 12611 ) for this particular case. 
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